function [F, ts, ws] = wft(f, window, support, NN, a, b, N, c, d, M)
% Window FT
% F(i,j) = F(ti, wj)

F=[];
ts = linspace(a, b, N+1);
ws = linspace(c, d, M+1); ws(end)=[];
x = linspace(support(1), support(2), NN+1);
for t = ts
    Ft = ffft(f(x+t).*window(x), support(1), support(2), NN, c, d, M);
    F = [F; Ft];
end
[ws, ts] = meshgrid(ws, ts);
end

